A dark matter interpretation for the ARCADE excess?B 
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The ARCADE 2 Collaboration has recently measured an isotropic radio emission which is signif- 
icantly brighter than the expected contributions from known extra-galactic sources. The simplest 
explanation of such excess involves a "new" population of unresolved sources which become the 
most numerous at very low (observationally unreached) brightness. We investigate this scenario in 
terms of synchrotron radiation induced by WIMP annihilations or decays in extragalactic halos. 
Intriguingly, for light-mass WIMPs with thermal annihilation cross-section, the level of expected 
radio emission matches the ARCADE observations. 

PACS numbers: 95.35.+d,95.30.Cq,95.85.Bh 



The detection of a non-gravitational signal of Dark 
Matter (DM) would be one of the greatest pillar of mod- 
ern physics, simultaneously confirming our views of cos- 
mology, astrophysics and particle physics. This possibil- 
ity might be not far ahead if DM is in the form of Weakly 
Interacting Massive Particles (WIMPs), which currently 
are the most investigated class of DM candidates [l[ . Sig- 
natures of this scenario include a multi-wavelength spec- 
trum associated to radiative emissions involving electrons 
and positrons generated in WIMP annihilations or decays 
(for a recent review on this topic, see e.g., Ref. Q). 

Recently, the balloon-borne experiment ARCADE 2 
(Absolute Radiometer for Cosmology, Astrophysics and 
Diffuse Emission) Q reported radio measurements of the 
sky temperature at frequencies ranging from 3 to 90 GHz. 
Observations have been performed on a region which is 
roughly an annulus centered at (1, b) = (70, 0) with radius 
and widths of 30 and 20 degrees, respectively 

An isotropic component can be isolated from the AR- 
CADE data by subtracting foreground Galactic emis- 
sion Q. Surprisingly, the level of the remaining flux 
(which has been interpreted in terms of extra-galactic 
sky temperature) is about 5-6 times larger than the to- 
tal contribution from the extra-galactic radio sources de- 
tected in current surveys 0] ■ Even extrapolating the 
source number counts to lower (unreached) brightness, 
such excess still remains. 

Most sources of systematic effects which could explain 
the ARCADE excess have been ruled out Q. An as- 
trophysical galactic origin appears to be rather unlikely 
(see discussions in Refs. 3 and Q). Indeed, free-free 
emission has been excluded based on the spectral shape, 
and diffuse Galactic synchrotron foreground is estimated 
using two different methods (namely, a co-secant depen- 



dence on Galactic latitude and the correlation between 
radio and atomic line emissions) , which agree well among 
each other. 

The observed isotropic temperature can be fitted by 
the CMB blackbody contribution plus a power law: 



T(v) = T Q 
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where T Q = 2.729 ± 0.004 K i is the CMB thermody- 
namic temperature. Performing analogous analyses on 
past surveys at 22, 45, 408, and 1420 MHz similar re- 
sults are obtained, and fitting all data simultaneously 
the ARCADE Collaboration derived a = -2.62 ± 0.04 
and T s = 1.19 ±0.14 K Q. 

Such level of cosmic radio background does not have 
an immediate explanation in standard astrophysical sce- 
narios. In Ref. [8j, radio supernovae, radio quiet quasars 
and diffuse emission from intcrgalactic medium and clus- 
ters (as well as a missed flux from well-known sources) 
have been considered, concluding that none of them can 
significantly contribute. A new population of numerous 
and faint radio sources (able to dominate source counts 
around /iJy flux) has to be introduced 0, Q . Ordinary 
star-forming galaxies with a radio to far-infrared flux 
ratio which increases significantly with redshift can in 
principle offer a solution. On the other hand, this possi- 
bility is strongly constrained by multi-wavelength obser- 
vations. Indeed, the radio to far-infrared emission has to 
be increased by a factor of 5 above what is observed in 
local galaxies , while current measurements show very 



1 Radio synchrotron sources are associated to young stars and 



mild evolution, at least up to z ~ 2 — 3 ll| . An explana- 
tion of the ARCADE excess through radiative emission of 
secondary electrons in star-forming galaxies would over- 
produce the gamma-ray background from pion decays 
[l2j. The same is true also for primary electrons unless 
such putative galaxies have extremely low gas density 
(and, in turn, low ratio of primary electrons to pions) or 
extremely efficient proton escape. The picture that seems 
to emerge from ARCADE measurements 0, @| and sub- 
sequent interpretations @, S 11 , li| suggests the need for 
a population of numerous and faint synchrotron sources 
generated by primary electrons with a hard spectrum and 
with no or very faint correlated mechanisms at infrared 
and gamma-ray frequencies. 

In our current understanding of structure clustering, 
any luminous source is embedded in a DM halo, and 
therefore extragalactic DM halos can be seen as the most 
numerous source population. The flux induced by WIMP 
annihilations/decays is predicted to be very faint. It is as- 
sociated to primary electrons and positrons generated as 
final state of annihilation/decay, and WIMP models with 
large annihilation/decay branching ratios into leptons in- 
duce hard spectra of e + je~ with very faint gamma-ray 
counterpart (and, of course, no straightforward thermal 
emission). Therefore, WIMP sources represent an ideal 
candidate to fit the ARCADE excess and in this Letter 
we quantitatively investigate such possibility. 

Assuming a one-to-one relation between the mass M 
of extragalactic DM halos and the intrinsic luminosity £ 
of the source, the total isotropic intensity per solid angle 
at a given frequency v is given by (for a more detailed 
derivation of equations considered in the following, see, 
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FIG. 1: Extragalactic radio background as derived by AR- 
CADE [f| , together with three possible interpretations of the 
low-frequency (< 10 GHz) excess in terms of WIMP annihi- 
lations (blue, green, and orange curves, see text for details). 
The astrophysical source contribution estimated from number 
counts (red line), the CMB contribution (black-dotted line), 
and a best-fit power-law of the excess (black-dashed line) are 
also reported @]. 



where R v is the virial radius of the DM profile p and / is 
the fraction of halo mass in substructures (with a = 1 and 
a = 2 for decaying and annihilating DM, respectively). 
In Eqs. © and ©, dn/dM(M) and dn s /dM s (f, M S ,M) 
denote the mass function of the DM halo and of substruc- 
tures, respectively. For synchrotron emission: 



e.g. Refs. [lalJJl): 
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where z is the redshift, H is the Hubble rate, M c is the 
minimum mass of an emitting halo, and the luminosity C 
is function of the redhsifted energy E = E V {1 + z) with 
E v = hv. The luminosity function, including also the 
contribution of substructures within the DM halo, can 
be written as: 



power [15| . B denotes the magnetic field, and v is the 
frequency of emission (as opposed to frequency of ob- 
servation in Eq.©). The electron/positron equilibrium 
number density n e is obtained solving a transport equa- 
tion for e~/e + injected by DM with an energy spectrum 
set by dN e /dE e [1J|. The source terms of this equation 
for annihilating and decaying DM are: 
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therefore tightly correlated to the star formation rate of galaxies. 
Empirically, a nearly linear far IR— radio correlation is reasonably 
well established (at least for high-brightness and low redshift), 
see e.g. the review in Ref. [lfj . 



where M x is the mass of the DM particle, (a a v) is the 
non-relativistic annihilation cross section and the de- 
cay rate. 

For what concerns the 'astrophysical' parameters, we 
focus on two benchmark cases, which are fairly realistic 
(for a detailed discussion of impact of astrophysical un- 
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FIG. 2: X-ray and gamma-ray fluxes for the three benchmark 
WIMP cases shown in Fig. Q] The CHANDRA @ bound 
in the X-ray band and the COMPTEL @ and FERMI [13] 
extragalactic gamma-ray fluxes are shown. 



certainties, see Ref. [lj]). We adopt the halo mass func- 



tion dn/dM from Ref. 16j, recent N-body simulation 



results for concentration of halos [l7[ , a DM distribution 
inside halos following a NFW profile 18 1, and the mini- 
mum halo mass is set to M c = 1O 6 M (model A, in order 
to consider only objects for which we can guess a reason- 
ably large magnetic field), and to M c = 10 _6 Mq (model 
B). The contribution from substructures is modeled such 
that / = 10% of the total mass is in substructures and 
dn s /dhi(M s ) oc 1/M S , which leads to a boost in the sig- 
nal of b. /. su b — 7 (model A) and no boost in model B. 
Magnetic field is assumed to be constant in space and 
time with magnitude B = 10/iG (model A) and B = 2/xG 
(model B), and e + /e~ are assumed to radiate at the same 
place where they are injected. 2 The normalization of the 
emission roughly decreases by an order of magnitude go- 
ing from B = 10 fj,G to B = 1 /j,G, increases by two orders 
of magnitude going from M c = 1O 6 M to M c = 1O~ 6 M , 
scales linearly with &./. su b, and is mildly dependent on 
halo mass function, concentration and profile of DM. 3 
The excess spectrum reported by ARCADE and de- 



2 Notice that for the DM candidate of specific interest here (i.e. 
light or inducing a soft spectrum of e+/e~), electrons are mostly 
emitted at GeV energies, namely are injected at energies relevant 
for radiation at GHz frequencies. Therefore they don't travel 
significant distances before radiating, while electrons emitted at 
larger energies would take some time to cool down, and diffusion 
and escape time would become much more relevant. 

3 For example, the case of an isothermal cored profile leads to a 
reduction in the intensity by, roughly, a factor of 1.5. 



scribed in Eq. ([T} is rather hard, and requires a hard 
electron/positron spectrum dN e /dE e . This can be pro- 
duced by DM scenarios with a large branching ratio of an- 
nihilation/decay into leptons. For illustrative purposes, 
we chose the fi + — fi~ channel. To reproduce the absolute 
normalization of the excess with a 'thermal' annihilation 
rate (a a v) = 3 • 10 _26 cm 3 s _1 in our benchmark model 
A, we need a WIMP with M x = 10 GeV (M x = 25 GeV 
in model B), which, even though induces a slightly softer 
spectrum than the best-fit power-law, provides a rea- 
sonable agreement with the data, as shown in Fig. [T] 
For this benchmark case we have x 2 /dof = 26.9/13. 
The actual best-fit for model A is obtained with a mass 
M x ~ 30 GeV, (x 2 /dof = 14.3/13), but at the price 
of increasing the cross section by one order of magni- 
tude. The fact that light DM, in the 10 GeV mass range, 
can fairly well reproduce the ARCADE excess, without 
the need of unrcalistically large DM overdensities is par- 
ticularly interesting, especially in light of recent claims 
of signals compatible with a DM interpretation from di- 
rect detection experiments (DAMA 0, CoGeNT 0, 



and CRESST [21]), that can be in fact accommodated 
with a - 10 GeV WIMP 0. In the case of the AR- 
CADE excess, the best option to explain the effect in 
terms of DM annihilation requires a light DM particle 
which annihilates mainly into leptons, and therefore that 
does not couple dominantly to quarks (coupling relevant 
to the direct detection scattering cross-section). Never- 
theless, it is not very difficult, from the model-building 
point of view, to foresee a model where a DM candidate, 
which annihilates mainly into leptons, still has relatively 
large scattering cross-section off nuclei. For a concrete 
example, see e.g., Ref. [23j. Note also that the radio 



emission in the Milky- Way halo induced by WIMPs fit- 
ting ARCADE data can either easily satisfy constraints 
(for cored galactic DM profiles) or be close to a possible 
detection in the central region of our Galaxy (for cuspy 
profiles) [24 1 . 

Similar conclusions on the viability of a light 'lep- 
tophilic' DM particle in explaining the ARCADE data 
can be also drawn in the decaying case: for a DM mass 
of M x = 10 GeV, the excess is reproduced if the lifetime 
is Td = 3 • 10 27 s (with a curve similar to the one shown 
for the annihilating case). 

Since the ARCADE excess can be explained by DM 
annihilation/decay in terms of sizable production of elec- 
trons and positrons, emissions of X-rays and gamma- 
rays by means of inverse-Compton processes on inter- 
stellar radiation fields (here we include CMB only) and 
direct production of gamma-rays from the DM particle 
annihilation (cither by production of neutral pions or by 



Final-Statc-Radiation (FSR)) arc present and have to 
be checked against available bounds. For the benchmark 
cases considered, these multi-wavelength constraints are 
easily satisfied, as shown for X- and 7-rays in Fig. [2] 

In Fig. [U we show also the case of a more 'classic' 
WIMP candidate with 100 GcV mass and hadronic anni- 
hilation channel (bb pair in the shown benchmark). This 
scenario is less appealing than previous cases. The spec- 
trum is relatively too soft in order to reproduce well the 
ARCADE data and the fit is worse (x 2 / do1 = 49.5/13) 
than for DM annihilating into leptons; note that the ex- 
cess is sizable up to at least 3 GHz (although not clearly 
visible in Fig. [T]duc the smallness of the error bars and 
the scale of the plot), so a viable explanation has to 
roughly overlap to the dashed best-fit curve up to those 
frequencies. Moreover, since now the DM mass is larger, 
the required boost factor is accordingly larger (by a fac- 
tor of 20 in this specific case), which can stem from a 
larger annihilation cross section related to a non stan- 
dard formation of DM relic density, from a Sommerfcrld 
enhancement, or from a larger contribution of substruc- 
tures with respect to what considered here (see e.g., [28 1 



for further details on possible boost factors). Heavier 
DM with hadronic annihilation/decay final states is also 
more strongly constrained by the 7-ray channel, as can 
be seen in Fig. O 

As a further analysis on the radio emission arising from 
light DM annihilation/decay, able to adapt to the AR- 
CADE excess, we show in Fig. [3] the differential number 
counts of sources at 1.4 GHz. If wc assume all substruc- 
tures to be unresolved, they mainly boost the signal of 
large and bright halos (since the latter host more subha- 
los). On the contrary, if all substructures are assumed to 
be resolved, counts drop much more slowly at low bright- 
ness. To bracket uncertainties related to the possibility 
of resolving substructures in the future, these two ex- 
treme cases are shown in Fig. [3] as solid lines, for the 
same 10 GcV DM benchmark of Fig. [T] As discussed 
above, the key point for our analysis is that in both 
scenarios the number of DM sources definitely becomes 
dominant over astrophysical contributions (AGN, star- 
forming galaxies) at the sub-/iJy level. The contribution 
of star-forming galaxies, which is dominant over AGN 
emission at low fluxes, decreases more rapidly (assuming 
FIR/radio correlation holds at all redshift), than the ex- 
pected contribution from DM, both in the case of resolved 
and unresolved substructures. From Fig. [3] we notice that 
the flattening at low brightness in current data, although 
it can be easily accounted for by standard astrophysical 
populations (for a review, see e.g. Ref. (lo||). neverthe- 
less could be well fitted by a DM model between the two 
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FIG. 3: Differential number counts for AGNs (dashed line), 
star-forming galaxies (dotted line), and the same 10 GeV 
benchmark DM model (solid lines) shown in Figs, [l] and [2] 
For DM, we consider a case such that all substructures are 
resolved, and an opposite case where all substructures are 
unresolved. For data and astrophysical models, see [l^] and 
references therein. 



extreme cases presented in Fig. [31 

In summary, we discussed the possibility that syn- 
chrotron emission induced by WIMP annihilations can 
account for the isotropic radio component measured by 
the ARCADE 2 Collaboration. Although galactic or 
extragalactic astrophysical interpretations of the excess 
cannot be excluded, they currently present some puz- 
zling issues @, 0, Q • Under reasonable assumptions for 
clustering, we found that light-mass WIMPs producing 
hard-spectrum electrons and positrons (as in the case 
of leptonic annihilation channels) in extragalactic halos 
with a 'thermal' annihilation rate can fit the excess and 
satisfy constraints at other wavelengths. A population of 
sources which can generally explain ARCADE measure- 
ments has to become the most numerous at brightness 
aroundujy, so it will be certainly studied in details by 
SKA [29(, and possibly also by its precursors, ASKAP 
El and MeerKAT " 



31] 



If the excess is due to extra- 
galactic DM, a clear discovery of a non-gravitational sig- 
nal of DM might be not far ahead. A dedicated study 
of closest and brightest (in terms of DM-induced signal) 
objects with current radio telescope (e.g., ATCA |32| and 
EVLA |33J) can start to probe this scenario in the near 
future. 
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